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ABSTRACT 

A recent survey of the Galaxy and M31 reveals that more than 90% of dwarf 
galaxies within 270 kpc of their host galaxy are deficient in Hi gas. At such 
an extreme radius, the coronal halo gas is an order of magnitude too low to 
remove Hi gas through ram-pressure stripping for any reasonable orbit distri- 
bution. However, all dwarfs are known to have an ancient stellar population 
(> 10 Gyr) from early epochs of vigorous star formation which, through heating 
of Hi, could allow the hot halo to remove this gas. Our model looks at the evo- 
lution of these dwarf galaxies analytically as the host-galaxy dark matter halo 
and coronal halo gas builds up over cosmic time. The dwarf galaxies — treated as 
spherically symmetric, smooth distributions of dark matter and gas — experience 
early star formation, which sufficiently heats the gas allowing it to be removed 
easily through tidal stripping by the host galaxy, or ram-pressure stripping by a 
tenuous hot halo (uh = 3 x 10~^ cm~^ at 50 kpc). This model of evolution is able 
to explain the observed radial distribution of gas-deficient and gas-rich dwarfs 
around the Galaxy and M31 if the dwarfs fell in at high redshifts {z ~ 3-10). 



1. Introduction 



In recent years, the Sloan Digital Sky Survey ( York et al. 200ol) has more than do u 



bled the number of known dwarf galaxies that orbit the Galaxv (IWillman et al. 



Zucker et al. I boOGbl : iBelokurov et al.ll2006l : IZucker et al. I l2006al : IBelokurov et al.l 



20071 : 



2005a b; 



Irwin et al. 



-2- 



20071 : IWalsh et all 120071 : iBelokurov et all I2OIOI : IWillman et all 120101 ) . These newly discov- 
ered dwarf galaxies, being among the most dark-matter dominated structures in the universe, 
help address the perceived imbalance between the number of predicted low m ass dark mat- 



ter substructure and tho se observed today, the "Missing Satellites Problem" (IKlypin et al. 



I999I : iMoore et al.lll999l ). As low mass substructures are the building blocks upon which 



galaxies are formed under the ACDM paradigm, the ability to probe dw arf galaxies provides 



invaluable knowledge into the history and formation of the Local Group (ITolstoy et al.ll2009 
Karlsson et al.ll201l[ ). 



Although the dark matter component can be tracked through the use of high reso lu- 
tion ALbody simulations Jpiemand et al.ll2007l . boosi ISpringel et aP boosi : IWetzell [201o[ ) . a 
complete model of the evolution of t hese dwarfs is still elusive. All of the known dwarfs 
show signs of old stellar populations (ITolstoy et al.ll2009[ ) with more of the star formation 
occurring in disc rete bursty periods of star formation — of order 25% of total star formation 
(ILee et al.ll2009l ) — compared to more massive galaxies which experience roughly continuous 
star formation. This bursty behaviour, with bursts s eparated by gigaye ars, does not neces- 
sarily require interactions to trigger star formation ( iBrosch et al.ll2004l ) with the blow-out 



and subsequent infall of neutral gas suff icient to create bursts in isolated dwarfs (IValcke et al. 
20081 : buillen fc Bland-Hawthornlboosl ). 



In a recent development, semi-analytic models of galaxy formation, combined with re- 
sults from high-resolution dark matter simulati ons, have been used to model the physical 
properties of a Galactic dwarf galaxy system ( iLi et al.l I2OIOI ). These models however do 
not yet a ccount for the underabundance of Hi detected i n these dwarf galaxies over the last 
40 years (lEinasto et al.lll974l : iGrcevich fc PutmanI l2009l . hereafter GP09). The dependence 
of the Hi deficiency on galactocentric radius is evidence of significant tidal and/or ram- 
pressure stripping of the dwarf galaxies which GP09 attribute towards close-in, potentially 
highly eccentric orbits, allowing the Hi to be removed by a hot halo surrounding the Galaxy. 
This removal for close pe ricentres (5, 50 kp c) has been explained by a combination of tidal 
and ram pressure forces (IMayer et al.ll2006l ). 



But we now show that unassisted stripping fails by an order of magnitude to explain the 
phenomenon of dwarf galaxy depletion, on scales of ~ 250 kpc as observed today, for any 
reasonable orbit families (©. In ^ we include the effects of early star formation and find 
that feedback-assisted stripping is essential to explain the observed effect. Using models of 
dwarf infall, we then compare results from our model to the observed properties of the Hi 
abundance in dwarfs today in ^ Finally, a discussion of the limitations of this model and 
what impact that will have on the results is given in ^ 
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2. Dwarf Galaxy Model 



A dwarf galaxy with an eccentric orbit will experience a wide range of extremes as 
it travels along its orbit, from a low velocity through sparse gas at its apogalacticon to a 
large velocity through dense gas with a strong external radiation field at its perigalacticon. 
All these environments must be modelled. We describe this wide range of environments by 
considering four scenarios: the orbital path of the dwarf around the host galaxy, discussed in 
§2.1( the heating and cooling of the gas present, discussed in §2.2( how tidal or ram pressure 
stripping removes gas from the dwarf in §2.3j and how this stripping can be infiuenced 
by stellar formation and feedback, detailed in §2.41 All four, to varying degrees, will be 
infiuenced by how much gas the dwarf galaxy has at the beginning of its orbit and the 
properties of the dark matter halo that envelopes the dwarf. 

The profile of the dwarf galaxy dark-matter halo is not a simple choice, with no solid 
consensus on which profile most accurately describes the low-mass subhalos present in dwarf 
galaxies. Observationally, the dark-matter profile of dwarfs ca nnot be distinguishe d between 
a cuspy profile and that of a cored, constant-density profile (IWalker et al.l 120091 ). Both of 
these profiles also arise through simulations w ith the cuspy Einast o profile arising through 
large A^-body, dark-matter only, simulations ( ISpringel et al.l 120081 ). With the addition of 
baryons, it has been argued that an even steeper central core — through the process of bary- 
onic contraction — ^can b e present than predicted by the dissipation less A^-body simulations 
(IBlumenthal et al.l Il986l : iDel Popolol l2009l : iNapolitano et al.ll2010l ). However, by coupling 
the baryons with the dark matter, and allowing dynamical and angular momentum trans- 
fer to occur, constant density cores with Burkert like profi les (iBurkertI Il995l ) are able to 
develop in simulation (|E1-Zant et al.ll200ll : iDel Popolol l2009l ). Since our major focus is on 
dark matter-halos with low gas fraction, whereby both baryon-dark-matter coupling and 
baryonic contraction will have a smaller impact, we ignore these competing factors and use 
the Einasto profile. 

The choice of profile for the host galaxy is less contentious, and is chosen to have an 
evolving Einasto profile. Although the host contains more baryons than the dwarf galaxies 
we are interested in, due to its much larger mass baryonic contraction processes have a lower 
net effect. Thus the evolution of the host galaxy is dominated by the infall of dwarfs onto 
the galaxy. The evolution of the host of this galaxy is assumed to be similar to that of the 
Milky Way, with a growth rate described by iMcBride et al.l (120091 ) leading to a final virial 
mass of 1.37 x 10^^ M (t, (for an Einasto m odel with the same circular velocity at virial radius 
as the NEW model in lSmith et al.l 120071 ) . to give a mass at redshift z of 



M{z) 



1.69 X 10"^^ / (1 + 1.172)Vfio(l 



Da— d2 + 0.96 

dz 



-7.87 



10^^ Up 



0! 



(1) 
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We assume that any cold gas is of a constant density, and gas in the warm and hot 
phases in the dwarf is in hydrostatic equihbrium with the density of the gas then given by 



n{r) = n(0) exp[—V{r)/c 



(2) 



where n(0) is the c entral density, is the speed of sound in the gas and V (r) is the Einasto 
potential given by ( ICardone et al.ll2005l : iNichols fc Bland-Hawthornll2009l ) 



V{x) = -f,(t)28-^/°exp(2/a)a-i+='/" [2^/"a-^/"7(2/a, 2x7a) - 7(3/a, 2x°/a)/x - l] , 

^ (3) 
with X = r/vsit) the scaled radius, rs(t) the halo scale radius, Vs(t) the characteristic scale 
velocity of the dark matter halo, a the Einasto free parameter, and 7(0, x) the lower incom- 
plete gamma function. 

The central density of the gas is calculated by assuming the calculated central density 
of Leo T in GP09 of nn = 0.45 cm~^ is comprised of gas in pressure equilibrium with a mass 
ratio 1 : 3 of cold gas to warm gas, matching the total ratio of masses within Leo T from 
which this value was calculated ( iRyan- Weber et al.ll2008l ). This gives a warm gas density of 
0.05 cm^^ and cold gas density of L5 cm~^. The hot phase central density is assumed to be 
in pressure equilibrium with the cold and warm phases. 



2.1. Orbit Of The Dwarf 



The orbit of the dwarf galaxy within the larger host galaxy is given by the shape of the 
potential. As the host galaxy evol ves, the halo scale radius and scale velocity will change 
according to ( ISternberg et al.ll2002l ) 



Vs{t) 



1/3 



1/3 



47rA(t)p„(t) 
47r 

V3G3A(t)p„(t) 



■^vir (^) 



1/6 



M. 



1/3 



/m,('Cvir(^)) 



(4) 
(5) 



where A is the characteristic overdensity of dark matter halos, pu is the mean baryonic 
density of the uni verse, Xv\r = Ryh/r,., is the concen tration parameter and fm is the mass 
profile as given in iNichols fc Bland-HawthornI ( 120091 ) . 



The characteristic overdensity is given by ( iBryan fc Norman! Il998l ) 

ISvr^ + 82{n{z) - 1) - 39{n{z) - 1)^ 



A 



(6) 
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and the mean baryonic density by 



Pu 



3H{z) 
SttG 



(7) 



For an Einas to profile, the concentration parameter, x^ir, also referred to as c, is 



flDuffv et al.ll2008r ) 



Xvir = 8.82(1 + z) 



0.87 



2 X 1012 h M. 



& 



-0.106 



(8) 



In this growing potential there is no angular dependence and therefore the angular 
momentum is conserved; such that 



3 ■ exp(2/a)t;2(t)«-i+3/" J^^j^^ 7(3/0, 2x'~''/a)dx' 



l/rl - l/rl 



(9) 



where is the per unit mass angular momentum squared and Tp and are the peri and 
apogalacticon radii that would occur if the potential was static. 

We present our results in terms of both orbit circularity and eccentricity in recognition 
of the fact that both are in common use when discussing CDM model. We define the 
eccentricity of the orbit in terms of the ratio of the apogalacticon to the perigalacticon, 
e = {ra—Tp) / [ra+rp], and the corresponding circularity to be 77 = \/l — e^. These parameters 
are not conserved throughout the orbit due to the growing potential, and when referred to 
we use the final state of these parameters. 

We also add dynamical friction as a non-conservative force, wit h a naagnitu de of the 
instantaneous dynamical friction timescale multiplied by the velocity ( iZhaol 12004 ) 



'-fric 



47rG'VMdw 



2.5 



4/3 + {v/v,^,f 



(10) 



giving the equations of motion as 



dV{r,t) AG'Kr'^puiz) 



dr 



— rt 



fric' 



e 



re 9 1 



(11) 

(12) 



with the additional term in the radial direction arising from the background contribution to 
the density of the dark matter, which becomes non-negligible beyond r ^ ^r^^^. 
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2.2. Heating And Cooling 

Throughout the passage of the orbit, the gas present within a dwarf galaxy will be 
exposed to a variety of radiation fields, with the most powerful being the radiation field from 
the host galaxy, the extragalactic background light and the stellar radiation field produced 
from high mass, short lived stars formed within the dwarf galaxy. To examine if the changing 
orbit itself is enough to produce the observed galactocentric effects, we first examine the 
system without the internal stellar radiation field. 

Without internal star formation, two contributions to the radiation field will dominate 
the heating and ionization of the dwarf galaxy, namely the extragalactic background field 
and the radiation field from the host galaxy which will dominate when the dwarf is within 
~ 100 kpc of the host gal axy's centre. For the e x traga lactic background field we use the 



redshift-dependent fields of iFaucher-Giguere et al.l (120091 ) and assume that all the radiation 



enters the dwarf galaxy radially. To a good approximation, the radiation field of the host 
galaxy is assumed to be a delta function at 13.6 eV wit h the number of photons impacting 



the dwarf galaxy at a distance r from the disk given by (IBland-Hawthorn fc Maloneylll999l ) 



, , 0.01 X 2.6 X 10^3 gpj^^^) 

^^^'") = 7^ SFR(0)' ^''^ 

with the 0.01 factor arising from a 6% escape fraction vertically from th e disk and isotropized 



over the halo. SFR(z) is the star formation rate at redshift z, given by (jJust &: Jahreifill2010[ ) 



qFPm-2 ^^ + ^-^^ Gyr)(11.7Gyr)3 

^^^^'^-^ — (t^ + [7.8 GjT]^y — ' ^^^^ 

where t is the time since disk formation, which we assume formed 12 Gyr before the present 
day. 

These two radiation fields combine to heat and ionize the cold gas that is present in the 
dwarf galaxy. The gas is ionized according to the absorbed radiation field, for which a radial 
radiation field yields 



e = ^ / J(r,z/,^)(l-exp[-iVH:a(z/)])dlogz/, (15) 

'log(13.6 cV/[27rrj]) 



h 



with i?dw the radial size of the cold gas in the dwarf and J(r, u, z) is the combined radiation 
field of the Galaxy and the extragalactic UV field at a distance r, frequency v and redshift 
z in ergs s~^ cm~^ Hz~^. 

These radiation fields will also heat the gas, with the heating (considering only H and 
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He) given by (IWolfire et al.lll995l ) 

flog(7.254xl0i^) 



XR 



h 



log(13.6 eV/[27rri]) 



J {r, u, z) {1 - exp[-Nui(r{u)])Eh{E,n J Tin) d log u, (16) 



where is the heating per primary electron given by IWolfire et aLl ( 119951 ) . 

The coohng, Af T, Z), is calcu l ated b y metal line cooling at a met allicity of O.lZg^ with 
the fit s provided by ISchure et al.l (120091 ) and He collision cooling by iDalgarno fc McCray 
Jl972h . 

As more particles are ionized than can be warmed, any excess ionizations (after recom- 
bination) are converted to heating with all 13.6 eV going to heat the gas. The transfer from 
cold neutral gas to warm ionized gas, is then given by 

TxR 



MxR 

Mcool 
^net 



8 







M, 



cold 



Meool - MxR - 



H,0,iu 

2 X 10-i°T^ 

^net - Mcool + MxR 







XR 



1.5kb{T^ - Tc)mH ' ^^'^^ 
where and are the temperatures of the cold and warm phase respectively, rin o^^ and 
Xedge,w IS the Central density and scaled radii of the edge of the warm phase respectively, and 
/gas is the Einasto gas distribution function. Correspondingly M„arm = — ^coid- 

Without any other source of heating (see §2.4p . the cooling prevents any significant 
heating in the model, keeping most gas as cold gas with only a small amount of warm gas 
in the centre. 



2.3. Tidal and Ram Pressure Stripping 



Gas in the warm or cold phase may be stripped through a combination of tidal forces 
and that of ram pressure arising from the orbit through the hot halo of the host galaxy. 

The tidal forces will not only potentially remove gas from the dwarf galaxy but will 
strip down the dark matter from the halo that surrounds the dwarf. T he tidal extent of the 
dwarf's halo, neglecting the contribution of baryonic mass, is given by (iHayashi et al.ll2003l ) 



m{Rt 



2 - 



r dM' 
M{r) dr _ 



M(r) 



'1^ 
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where Rt is the tidal radius of the dwarf, and r is the orbital radius of the dwarf from the 
host galaxies centre. For an Einasto profile, this tends not to have an analytic solution, and 
for reasons of speed is only solved initially, before using the derivative to calculate the change 
in tidal radius over time. Any gas that extends above the tidal radius is considered to be no 
longer bound, and hence lost, in order to simplify calculation. A more extensive treatment 
involving re-accretion onto dwarfs requires hydrodynamical treatments that are now under 
way. Within our models, tidal stripping does not produce significant gas loss with Galactic 
and extragalactic ionizing fields keeping most gas deep within the tidal radius. 

Working at the s ame t ime as the tidal stripping is ram pressure stripping. Here we use 
the [McCarthy et al.l ( 120081 ) formalism whereby the ram pressure removes shells of gas from 
the dwarf galaxy, assuming that the gas distribution behaves as an isothermal sphere beyond 
the cloud edge if 



Pgai(r)f > - 



R 



(19) 



This formalism, originally used for a hot halo, can be carried over to the warm and cold 
phases by assuming spherical symmetry and a p oc density drop off outside the cloud 
edge. 

We note that when u sing the realist i c appr oximation a ~ Vchc /V2, this is Svr times more 
difficult to strip than the iGunn &: GottI (119721 ) equation used by GP09, i.e. 



an. 



■gas 



(20) 



If the stripping does occur, then the rate of stripp ing will be dependent upon the speed 
of the shock induced in the gas (IMori fc Burkertll2000l ). such that 



Vfs 




(21) 



Here we assume that the halo stays spherically symmetric at all times, and hence the rate 
of the change in radius will be half the speed of the shock, i.e. R = 0.5vfs- 

The density of the host galaxy halo will be dependent upon the temperature of the halo. 
We use an isothermal halo with the temperature set at the virial temperature of the halo at 
any given redshift. We set an unchanging central density to give a density of 3 x 10~^ cm~'^ 
at 50 kpc at the present da y. These values are consistent with isothermal hal o parameters 
from iBattaglia et al.l (120051 ) . the values used by iBland-Hawthorn et al.l (120071 ) (which uses 



T = 1.75 X 10^ K, nn = 2 X 10 ^ cm ^ at 55 kpc), and the high-entropy Galactic halo model 



of iKaufmann et all (120091 ). 
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For this stripping, the phases are assumed to be independent, and the cold phase can 
be removed at the same time as the warm phase if the ram pressure is large enough. The 
rate of change of mass is then simply 



M = 47rrVgas(i?)A 



(22) 



With only Galactic and extragalactic ionizing fields, ram-pressure, like tidal stripping, 
does not produce large changes in gas mass even over cosmic time. Even with more favourable 
conditions to ram pressure stripping, such as the Gunn-Gott criterion and a lower tempera- 
ture halo (resulting in a denser inner portion), only a small amount of gas is lost — less than 
10^ M0 — corresponding to the cold gas heated by the Galactic and extragalactic UV fields 
to a warm phase and then lost to ram-pressure stripping. 

These results are in conflict with much more substantial gas loss has been measured 
in hydrodynamical simulations wit h ram pressure alo ne able to strip the gas from close in 



pen 



< 5 kpc), low rnass h alos (IMayer et al.l 120061 ). However, at larger perigalacticons 



dwarfs in iMayer et al.l (120061 ) are still able to retain much of their centrally concentrated 
gas, suggesting another form of heating is required to allow the removal of gas from dwarfs 
with larg er perigalacticon s (e.g. the perigalacticon of Sextans and Draco is expected to be 
> 50 kpc lLux et aDboioh . 



2.4. Star Formation and Feedback 



A dwarf consisting of cold, dense gas is resistant to tidal and ram-pressure stripping 
with only a thin skin, ionized by the Galactic of extragalactic UV fields being removed. 

Early star formation however, will heat this cold gas, raising it in the potential well and 
making it more easily stripped. This warm ionized gas, at the same pressure as cold gas will 
occupy ~ 60 times the volume, and being much less dense will become even easier to strip 
than its height in the potential well would indicate. 

The star formation in dwarfs is considered to consist of periods of low-level star forma- 
tion, during which short bursts are induc ed which i n crease the star formation by a factor 
of 3, consistent with dwarfs observed by iLee et al.l (120091 ). We assume these bursts are 
triggered by perigal acticon passages, induced by shocks created through tidal interactions 
(jPasetto et al.ll2010l ). There i s also evidence tha t bursts may be triggered by the re-accretion 
of heated and expanded gas (IValcke et al.ll2008l ). however, much of this gas will be stripped 
while in the warm phase and far from the centre preventing the re-accretion and subsequent 
starburst from occurring. 
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The base star formation rate is taken to be similar to the dwarfs that surround M31, 



with a star formation rate (IKaisin &: KarachentsevI l2006l ) 



log(SFR Mo yr-^) = 1.4 log Mm + k, (23) 

where k is a constant calculated by assuming that the gas will be completely depleted at a 
time t = 5/Hq 70 Gyr. For an initial gas mass of 5 x 10^ M©, k = —13.62. 

The UV-X-ray spectrum produced by the star formation in the d warfs — which is respon - 



sible for the majority of heating — was calculated with Starburst99 (iLeitherer et al.lll999l ): 
when a burst occured the change in the spectrum was calculated by summing many short 
bursts together to produce an approximately continuous change in the spectrum. Above the 
Lyman limit, this spectrum (in ergs s~^ cm~^ Hz~^ sr~^) and its changes during a burst 
was well fitted by 

log J, = 2 + Ji + log[SFR Mq yr-i] - 3.54 log z/, (24) 

where Jj changes linearly over a 5 Myr period at the beginning and end of each burst from 
Ji = 96.15 (ergs cm~^ Hz~^ sr~^) at the low level continuous star formation rate to 

Ji = 96.63 (ergs s~^ cm~^ Hz~^ sr~^) during a burst. 

The addition of star formation will also result in supernova superheating the gas well 
above the temperature of the warm phase Any power from the supernova is assumed to go 
into heating gas into a hot (10^ K) phase. In reality, the energy from the supernova will also 
go into expanding the cold and warm phases, making them more susceptible to stripping, 
but due to the hmitations of a smooth density profile we assume this energy is instead used 
to heat the gas. 

We find the power output of a supernova for continuous, non-bursty star formation is 
given by 

log[PsN erg s"^] = log[SFR M© yr'^] + 40.0. (25) 

After a burst has commenced, the power output by the supernovae increases over 38 Myr, 
with the power over this period of increasing output given by 

log[PsN erg s~^] = log[SFR Mq yr~^] + 2 + 25.6 tan-i[24(0.1t + 1.7)], (26) 

where t is the time in Myr since the burst commenced. 

During the increased star fo rmation associated with the bursty period (lasting < 10^ yr 



Sharp fc Bland-Hawthornll2010l ) the power output is 



log[PsN erg s"i] = log[SFR Mq yr'^] + 42.0. (27) 
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After the burst has ended and star formation drops back down to the continuous star 
formation rate, the power output subsides as massive stars die off. Over a period of 40 Myr, 
the power output dechnes according to 

log[PsN erg s-i] = log[SFR yr'^] + 42.0 - (0.0125tO, (28) 

where t' is the time since the starburst has ended in Myr. 

In this simphfied model, we assume that any power from the supernova is dissipated in 
the gas and goes into heating the cold and warm gas at 100% efficiency and split between 
the cold and warm phase in proportion to their mass. The mass loss rate will then be the 
power dissipated into that phase divided between the energy difference between that phase 
and the hot phase. 



1.5fc,[(T, - TJ + (T, - T^)^] 

Although the efficiency of supernova heating is typically much lower, there is only a small 
amount of material that is heated into the hot phase, with much more material being lost 
through the warm phase. In this limit, the snowplough phase of supernova expansion, where 
a large portion of energy goes, may have a noticeable contribution by raising gas outside of 
the potential well and we keep the efficiency at unity to partially compensate for this effect. 

Under these conditions the internal star formation greatly exceeds the extragalactic UV 
background at early times, making it the dominant source of heating. The cooling of gas 
from the warm phase to the cold phase can approach the level of UV-X-ray heating, however, 
the amount of warm gas that is able to be held on by the halo is limited, preventing most 
gas from recooling. The relative rates of heating by the combined radiation fields (internal 
star formation. Galactic and extragalactic) and the cooling rate is shown in Fig. [1] for 
a model with perigalacticon of 80 kpc and circularity of 0.9 (eccentricity of 0.44). Very 
little gas is moved to the hot phase from supernova bursts, much more gas is heated to the 
warm phase which quickly extends beyond the tidal radius, resulting in large amounts of gas 
loss from tidal stripping. Only in more massive dwarfs with correspondingly larger central 
densities will sufficient quantities of warm gas be able to cool back to a cold phase to allow 
for significant supernova heating. 
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Heating and Cooling in a dwarf galaxy 
Redshift 
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Fig. 1. — The relative rates of cooling for a dwarf galaxy finishing with perigalacticon of 
80 kpc and circularity of 0.9 (eccentricity 0.44). Initially the cooling is minimised by warm 
gas being forced outside the tidal radius of the dwarf, preventing it recooling onto the dwarf, 
after the first pericentre passage, the cooling and heating reach equilibrium, until the next 
pericentre, when the warm gas drops below 10 Mq and is no longer tracked. The drop in 
cooling around 7.5 Gyr is an artifact of the method of integration when heating and cooling 
are no longer in equilibrium and does not change any results. 
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2.5. Method and Initial Conditions 

The gas and dark matter masses were tracke d by solving; the ti me dependent equations 



equations flT|)- fl28|) — with the odeint routine (jPress et al.l Il992l ) and encoded in python. 
Once a gas phase dropped below 10 Mq, it was no longer tracked and assumed to no longer 
have any gas in that phase. 

The dwarfs were initiated with 5 x 10'' Mq for the cold phase and 100 Mq for the warm 
and hot phases. The dark matter halo of each dwarf galaxy was set with a virial mass of 
1 X 10^ Mq for models beginning at z = 1 and z = 3 and a virial mass of 3 x 10^ Mq for 
z = 10. For all starting points, the predict ed mass within 300 p c is about 10^ Mq, consistent 



with the mass inferred for dwarf galaxies fIStrigari et al.ll2008l ). 



The orbit of the dwarf galaxy was set so that it ended at apogalacticon, with the initial 
conditions found by integrating backwards without dynamical friction initially, and then 
with dynamical friction along the previous orbit until convergence of < 5% at all radii was 
found. The model was run for dwarf galaxies that at their present day orbits will have a 
pericentre of between 15 and 250 kpc, and a circularity of between 0.2 and 0.9. We assume 
each point is similar to those around it, and each model was at the centre of a 5 kpc, 0.05 
circularity box, for those with a pericentre Vp < 50 kpc, a 10 kpc, 0.1 circularity box for 
pericentre 50 < Vp < 100 kpc and a 25 kpc, 0.1 circularity box for 100 < < 250 kpc, 
for a total range of 12.5-262.5 kpc and 0.175-0.925 in circularity. A sample of these orbits 
(Fig. [2]) illustrates the effect of varying eccentricity and pericentre in a growing host halo 
with dynamical friction. 



3. Results 

Beginning at z = 1-3 allows most halos to retain a large portion of their cold gas, while 
beginning at higher redshift 2; = 10 means most halos have been stripped of nearly all their 
cold gas. In all cases halos lose over 90% of gas mass to warm phase stripping resulting from 
early heating by internal star formation. As could be expected, highly circular orbits at low 
pericentre show the strongest gas loss, with the final mass displayed in Fig. [31 This is large 
gas loss attributable to an on-average larger number of bursts and the increased halo density 
at low pericentres. We note that finishing at apogalacticon will always result in later bursts 
for the same pericentre, but the results remain if the models begin at apogalacticon in which 
case the bursts will always be earlier. 

A large portion of the parameter space exists in which a halo will always be within 
270 kpc and retain over 1 x 10^ Mq of gas. Numerous other halos will spend a portion 
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Fig. 2. — A small subset of orbits of dwarf galaxies of mass Myir = 3 x 10^ Mq beginning at 
z — 10. The effects of a growing potential is mostly easily observed with the early orbit of a 
low pericentre, high eccentricity dwarf (r^ = 25 kpc, e = 0.92), whose changing orbit cannot 
be accounted for by dynamical friction alone. 
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Fig. 3. — The final neutral hydrogen mass for dwarf galaxies beginning at z = 1. The 
figures show equivalent information: (a) is the orbit eccentricity vs. the radius of closest 
approach (perigalacticon) today; (b) is the orbit circularity vs. perigalacticon (see §2.ip . The 
rectangular box in both figures shows the range of parameter space explored in this work 
(see The dashed line corresponds to an apogalacticon radius = 270 kpc, the radius at 
which dwarf galaxies suffer strong Hi depletion in M31 and the Galaxy. The greyscale (color 
scale in online version) depicts the amount of Hi retained by a dwarf galaxy observed today. 
Although the gas in each phase was tracked down to 10 we restrict the range above to 
10^-10^ Mq for readability. Galaxies with < 270 kpc can only retain Mhi < 5 x 10^ Mq. 
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of their orbit within this 270 kpc hmit, but as the model finishes at apogalacticon, these 
subhalos may have experienced less pericentres than would be expected if they were within 
270 kpc today. 

For those subhalos which begin at z = 3 (Fig. H]) no halos are always within 270 kpc 
and maintain over 1 x 10^ M© of Hi; although some halos should be detectable with over 
1 X 10^ Mq of gas remaining that will always be within 270 kpc. 

The low circularity, high pericentre drop in gas mass is a result of the model finishing at 
apogalacticon. This effect is much more pronounced for halos beginning at 2; = 10 displayed 
in Fig. ISl here the initial velocities may exceed 1000 km s^^. The long time in orbit, and 
lower virial mass of the dwarf galaxies in this run means that nearly all halos that retain 
gas will not have an apogalacticon within 270 kpc, that is, any observed halos with gas that 
began at z = 10 will only be observed within 270 kpc for a portion of their orbit. 

In all models, at early redshifts the high rate of star formation produces large quantities 
of warm gas which is removed by tidal forces existing on the dwarf galaxy, and warm gas 
at the edge of the tidal radius is then quickly removed by stripping of the hot halo. This 
initial period of extremely high mass loss accounts for over 90% of all mass loss as seen in 
Fig. [61 This initial period of gas loss is followed by a steady state, seen more easily in the 
mass loss rate versus time in Fig. [TJ is reached where cold gas is not warmed rapidly, and a 
small amount of warm gas can survive stripping, with most of the gas being stripped from 
dwarfs which are approaching their perigalacticon and experiencing bursts of star formation 
and the increased impact of the Galactic UV field. 

At low redshifts, the rate of stripping plummets as the initial conditions require that 
dwarf galaxies be at their apogalacticon at z = 0. Distributing the dwarf galaxies in phase 
along their orbits today has a minimal effect on the results presented in Figs. |3H51 The steep 
rate of gas loss for dwarf galaxies beginning at z = 10 is also a byproduct of our requirement 
that dwarfs finish at apogalacticon today. This leads to the excluded region (marked 'RAD' 
in Fig. [5]). The reason is that the dwarf is forced to orbit at unrealistic velocities to get back 
to its apogalacticon within an evolving potential well. This artifact has only a minimal effect 
on our results. 



4. Comparison with observations and completeness 

We now determine the fraction of dwarf galaxies retaining gas within 270 kpc of the 
Milky Way and M31. In order to estimate this, the completeness of the parameter space 
searched needs to be calculated. We assume that the infalling halos follow a distribution 



-17- 



Cold HI at z=0 for dwarf galaxies beginning at z= 3 (Mp^ = 1x10^ M^, IVIgg^ = 5x10^ M^) 
Perigalacticon versus eccentricity Perigalacticon versus circularity 




50 100 150 200 250 50 100 150 200 250 



Perigalacticon (kpc) Perigalacticon (l<pc) 

Fig. 4. — The final neutral hydrogen mass for dwarf galaxies beginning at z = 3 (Color 
scale available online). See caption to Fig. [31 Galaxies with < 270 kpc can only retain 
Mhi < 10^ Mq. 
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Cold HI at z=0 for dwarf galaxies beginning at z= 10 (Mp^^ = 3x10^ M^, Mg^^ = 5x10^ M^) 
Perigalacticon versus eccentricity Perigalacticon versus circularity 
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Fig. 5. — The final neutral hydrogen mass for dwarf galaxies beginning at 2; = 10 (Color 
scale available online). See caption to Fig. [31 Galaxies with < 270 kpc can only retain 
Mhi < 10^ M©. 
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Fig. 6. — The total dwarf Hi loss vs cosmic time due to ram pressure stripping. (Tidal 
stripping from more massive galaxies, e.g. the Magellanic Stream, is likely to lead to higher 
mass loss rates at later epochs.) This is calculated by applying the weighting of each dwarf 
(see m]) and multiplying by the number of known dwarf galaxies orbiting the Galaxy (~ 33) 
today. We note that the amount of dwarf galaxies seen today is a small fraction of the 
total dwarf galaxies that have existed over time, with many destroyed by the Milky Way 
previously and others unobserved, the total mass loss above will similarly be only a small 
fraction of the total mass loss from dwarfs over time. The majority of this mass loss occurs 
very early on as the large amount of cold gas allows rapid star formation to occur, quickly 
heating and expanding the gas. This low density gas far from the centre of the potential 
well is easily removed quickly, before a plateau of mass loss is eventually reached. 
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Fig. 7. — The total dwarf Hi loss rate vs cosmic time due to ram pressure stripping. (Tidal 
stripping from more massive galaxies, e.g. the Magellanic Stream, is likely to lead to higher 
mass loss rates at later epochs.) This is calculated by applying the weighting of each dwarf 
(see Q and multiplying by the number of known dwarf galaxies orbiting the Galaxy (~ 33) 
today. We note again that the mass loss rate will only be a small fraction of the actual 
mass loss rate over all time as it does not include dwarfs that have been destroyed by the 
Milky Way or are presently unobservable. The period of greatest mass loss occurs when 
initial star formation heats the Hi allowing its removal from the dwarf galaxy before it cools 
down, a consequence of our model is that this occurs at the same time for all dwarf galaxies 
exacerbating this mass loss. This time period of greatest star formation rate is also when 
the super nova-driven winds are most importarit at r emoving gas but unlikely to explain Hi 
depletion (jMartirull998l : iBland-Hawthorn et al.ll201ll ) and are only accounted for in a simple 
way (see This figure assumes that the dwarfs outside our parameter space experience 
similar loss rates to those within. 



- 21 - 



identical to that of IWetzell ( 120101 ). which used a dissipationless A^-body simulation with 
ACDM cosmology (fi^ = 0.25, = 0.75, h = 0.72, n = 0.97 and o-g = 0.8) with a particle 
mass of 1.64 x 10^h~^ Mq. We note that although dwarf galaxy sized halos cannot be 
represented in this model at 2; = we extend the distribution on the basis that for higher 
masses there exists no dependence on mass. For a Milky Way sized halo at 2; = these 
distributions are 



df 

d(rperi/rv 



^ = 5.0r7^°^(l - vf'', 
drj 

3.6exp{-[3.4(rperi/rvir)]°-''} 



(30) 
(31) 



Taking the circularity and peri centre to be independent, the halos in the parameter 
space searched will comprise ~ 76% of all halos with a pericentre < 270 kpc and ~ 70% of 
all halos, with the distribution within the parameter space shown in Fig. [HI We split the 
distribution into four main regions: Regions A and B consists of those halos which were 
within our parameter space and always/sometimes within 270 kpc respectively. Regions 
C and D (halos with circularitys rj > 0.925) consists of those halos that were outside our 
parameter space and always/sometimes within 270 kpc respectively. Regions C and D were 
excluded from the parameter space as the large number of pericentres and low tidal gradients 
would give an unrealistic amount of bursts as well as the large number of pericentres being 
computationally expensive. The shaded region of low circularity/high eccentricity consists 
of halos which were outside our parameter space and spend so little time within 270 kpc 
that they were ignored. The shaded region of low pericentre is halos that would be stripped 
due to their very close passages, but the no disk assumption begins to break down as well 
as the increased computational cost means these halos were not modelled. 

In our model all dwarfs end at apogalacticon, however, the results are substantially the 
same if you begin them at apogalaction (for low redshifts where the orbit does not change 
significantly over time), we hence assume for analysis that they are at a random point along 
their orbit, and that this orbit is elliptical, as opposed to the standard Rosetta orbit an 
Einasto profile will produce. Under these assumptions, a dwarf galaxy with a pericentre of 
Tp < r' and eccentricity e = ^J\ — rf will be within r' kpc with probability 

1 / TT 

Pir < r ) = — a — arctan 

^ ^ TV \2 

a = ^ (r'/rp - 1)(1 - e)(l + e - (rVr,)(l - e)). 
Within each point in the model we then weight it by its contribution to all halos, with 



(r7r,)(l 



a 



(32) 
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Fig. 8. — The distribution of subhalos in pericentre-eccentricity (pericentre-circularity) space 
as a percentage of all subhalos with perigalacticon < 270 kpc at z = for a Milky Way 
sized galaxy, assu ming that circ ularity and pericentre are independent. Calculated using 
distributions from [Wetzel ( 120101 ). 
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the weight then 



w(r 



/ F[r < r ) ——ariar, 



(33) 



dr] dr 

where x is half the size of the box in kpc, and y is half the size of the box in circularity. 

We assume that all halos with a pericentre < 12.5 kpc cannot retain gas, and that 
subhalos with a circularity rj > 0.925 — halos from regions C and D — will have the same 
amount of gas, as halos with the same pericentre and a circularity of r/ = 0.9. Halos with a 
circularity 1] < 0.175 were assumed to contain gas, however, these halos have very eccentric 
orbits and hence will normally have a near zero weight. 

We then calculate what fraction of observed dwarf galaxies contain gas, looking at all 
galaxies below a given galactocentric radius in 100 kpc bins out to 2000 kpc and compare 
this to the fraction of galaxies — with confirmed detections and excluding the SMC and LMC 
discussed below — that contain gas from GP09. We use three cutoffs to distinguish between 
gas-deficient and gas-rich dwarfs. A lower cutoff of 10^ Mq is above the upper-boundary of 
a number of GP09 galaxies A mid-range cut-off of 10^ Mq is above the upper-boundary of 
even more galaxies in the GP09 sample. Finally an upper cut-off of 10^ M©, this cut-off is 
below the Hi mass of any confirmed gas-rich galaxy in the GP09 sample. These cut offs do 
not accurately represent the physical constraints on observations, with the observed being a 
combination of low cut-offs for nearby Galactic satellites to higher-cutoffs for far-off or M31 
satellites. 

Under our model, th e LMC and SMC being much more massive than t he typical dwarf — 
Mtm c: > 2 X 10^° Mq (ISchommer et al.lll992l ). MgMc > 3 x 10^ Mq jHarris fc Zaritskv 
20061 ) — will be able to retain significant quantities of warm gas against ram-pressure strip- 
ping. This protection from stripping would allow a significant quantity of gas to survive 
until the present day, with only tidal interactions — potentially those between the LMC and 
SMC — removing material from either dwarf galaxy. Although tidal stripping is included 
in our model, it is too simplistic to model the multi-body system present in the Magellanic 
clouds, with hydrodynamical simulations required to simulate the LMC-SMC-Galaxy system 
jLin et al.lEggsl : iBland-Hawthorn et al.lboOTt iNidever et aPboOsI : ISesla et al.lboiol ). 



For models beginning at z = 1 (Fig. |9]) we see a large over estimation for all gas cut-off 
masses. For models beginning at z = 3 (Fig. [TU]) again the predicted fraction exceeds the 
observed fraction of dwarfs, but predicts a lower fraction than models beginning at 2; = 1. 
Discovery of new gas-rich dwarfs far from the Galaxy may raise the observed fraction, but 
this may not explain the over- abundance predicted close in, where gas-rich dwarfs would be 
more detectable than further gas-deficient dwarfs due to more ongoing star formation. 
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Fraction of dwarf galaxies with gas today, Zg(3|^=1 
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Fig. 9. — The fraction of dwarf galaxies within a radial distance that contain potentially 
detectable amounts of Hi at z = ver sus radius for dwarf g alaxie s that began at 2; = 1. We 
note that the low radius spike in the iGrcevich fc PutmanI ( 120091 ) data (GP09 data) occurs 
due to a single galaxy — NGC205, a close in satellite of M31. Excluding this inner region, all 
three cutoffs greatly exceed the observed fraction of galaxies that retain gas. 
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Fig. 10. — Fraction of dwarf galaxies within a radial distance that contain potentially de- 
tectable amounts of Hi at z = versus radius for dwarf galaxies that began at z = 3. 
Excluding the inner region (see Fig. E) the cut-offs of 10^ Mq and 10^ exceed the ob- 
served fraction of galaxies that retain gas. The 10^ Mq cutoff closely follows the cutoff before 
exceeding it at large radius, this is likely a selection effect (see text). 
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For dwarf galaxies that begin their orbits at z = 10 (Fig. [TT 



the Mgas > 10^ M0 and 
Mgas > 10^ M0 cut-offs again overestimates the fraction of halos which will contain gas. The 
-^gas > 10^ Mq cut-off closely follows the shape of the GP09 sample, and is a lways below 
this fraction. At any distance we would expect the iGrcevich &: PutmanI (j2009[ ) fraction to 
be slightly higher than reality, as the gas poor dwarfs being difficult to see may not be fully 
accounted for. That the toy model represents the gas fraction best for halos beginning at 
2; = 10 a dds support to the idea that the Milky Way dwarfs were formed at or before this 
redshift JLux et al.ll201ol ). 



5. Summary &; Conclusion 



This simple toy model reproduces the observed fraction of gas-rich dwarfs, however, 
there are several factors that were assumed to be negligible which could affect the amount 
of gas that survives. The accretion of gas onto dwa rfs is unaccounted for in this model. 
The prospect of low redshift accretion (IRicottil 120091 ) in particular would greatly increase 
the chance of a dwarf surviving with gas to the present day, and if included for all dwarfs 
would likely result in an overestimation of the gas fraction even for dwarfs beg inning at late 
redshifts. The effects of dust through photoelectric heating and cooling is also ignored. Due 
to the low metallicity environments that dwarf galaxies typically have, this effect will be 
smaller than in larger galaxies, but may still be an important source of heating or cooling. 
We also assume that the early extragalactic UV field is uniform in space, a clumpy radiation 
field around the time of reionization, may have a large impact on the amount of gas that 
remains cold and protected from ram pressure stripping, this early gas loss could greatly 
impact the survival around the earliest pericentre passages. The use of a smooth medium 
for the gas — compared to a more realistic fractal medium — minimises the cooling of warm 
gas and allows it to extend beyond the tidal radius. Even with a large filling factor, the gas 
will cool much more quickly via metal line cooling than in the smooth medium minimising 
the large gas loss at the beginning. Potentially the biggest limitation is the assumption that 
all of the supernova energy goes into heating the gas to an extremely hot state, much of 
this energy likely goes into raising cold and warm gas out of the potential well of the dwarf, 
allowing it to be much more easily stripped. In particular this will predominantly impact 
the lowest gas masses, where the smooth gas assurnption versus a fractal medium is more 
likely to have a large impact ( Bland-Hawthorn et al.l 120071 ). 



Due to the large initial gas loss, our star formation rates do not represent that of the 
majority of dwarfs, with most stars forming early in the dwarfs life, as opposed to a roughly 
continuous star formation rate with some bursts (jWeisz et al.l |20 111 ). This makes our model 
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Fig. 11. — Fraction of dwarf galaxies within a radial distance that contain potentially de- 
tectable amounts of Hi at 2; = versus radius for dwarf galaxies that began at z = 10. Again 
the fraction of close-in dwarfs with gas observed exceeds that of the model (see Fig. [9]), at 
higher radii are bracketed by the Mhi = 10*^ and Mhi = 10^ M© cutoffs, discussed further 
in gl 
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more suited to explaining dwarfs with the se early periods of star formation where a majority 
of stars are formed, e.g. BK5N, KDG52 dWeisz et allboilj ). 



Even with these limitations, we believe that the toy model provides strong support 
towards internal heating due to early star formation, allowing the gas to more easily stripped 
due to the larger scale heights and lower density of warm and hot gas. In particular we were 
able to reproduce the fraction of dwarf galaxies th at retain Hi assu ming that the dwarf 
galaxies infall at a redshift of z = 10 consistent with iLux et al.l (120101 ). 



The covering fraction of warm gas stripped from the dwarfs is expected to be low, with 
the contrails containing approximately the same amount of war na gas as exists inside the 
cold streams that cover ~ 2% of the projected area of galaxies (iFaucher-Giguere &: Keres 



2010l ). This gas breaking away in small clumps of warm material will form a warm compo- 
nent of the hot halo whi ch is stabilised against heat conduction from the halo by cooling 
Jvieser fc Hensleil l2007bl lah 3efore falling onto the Galaxy as a warm rain. Even with the 
majority of dwarf galaxies entering at early redshifts, the large timescale required to com- 
plete large pericentre orbits (13 Gyr for Vp = 270 kpc) means that gas from any late infalling 
dwarfs — such as those from z = 3 or z = 1 — will only be falling onto the disk of the Galaxy 
today. A search for this gas in Ha around local galaxies is now being undertaken using the 
Maryl and- Magellan Tun able Filter ( iVeilleux et al.ll2010l ) and the Grantecan Osiris Tunable 
Filter (ICepa et al.ll2003l ) where the expected emission measure of the gas (< 0.1 cm~^ pc) 
may be just large enough to be detectable. 
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